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Abstract 



5h , 

Qh, The weak cosmic censorship hypothesis can be understood as a statement 

< 

, that there exists a global Cauchy evolution of a selfgravitating system outside 



an event horizon. The resulting Cauchy problem has a free null-like inner 
boundary. We study a selfgravitating spherically symmetric nonlinear scalar 



(N 
> 
in 
o . 

, field. We show the global existence of a spacetime with a null inner boundary 

© : that initially is located outside the Schwarzschild radius or, more generally, 

On 



outside an apparent horizon. The global existence of a patch of a spacetime 
that is exterior to an event horizon is obtained as a limiting case. 



1. Introduction 

The cosmic censorship hypothesis || can be informally stated as singularities are hidden 
inside black holes. Various attempts to formalize that statement have led to a collection of 
results that can be assembled into two categories, one of which may be called "geometric" 
and the other " dynamic" . 

The first approach, reviewed for instance in || and ||, seems to be strongly influenced 
by the post-Leray 0] notion of the global hyperbolicity. The latter requires (in addition 
to the standard hyperbolicity condition as stated for instance in ||) stringent smoothness 
properties of coefficients of a hyperbolic operator as well as a causality condition. In a 
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related stream of research a considerable effort is put into the examination of geometric 
quantities in various models, in order to define a notion of a singularity and then, to identify 
those singularities that must be enclosed by event horizons. A singularity can be under- 
stood as a point of space-time having divergent curvatures; a cosmic censorship hypothesis 
would require that all geodesies originating at a singularity remain inside a black hole. The 
post-Leray apparatus of ideas makes it plausible, however, to associate singularities with the 
existence of geodesies having a finite length |J . On the other hand, an increasing number of 
counterexamples to resulting versions of the cosmic censorship hypothesis has given inspira- 
tion to various notions of truly "singular" singularities, of which I mention a subclass of the 
so-called strong singularites |7j that is associated with the existence of apparent horizons 

The dynamic approach has been initiated, to my knowledge, by Eardley and Moncrief 
||. It bases on the notion of a Cauchy problem as formulated, for instance, in ||. A 
weak version of the cosmic censorship hypothesis can be formalized as follows, in the case 
of asymptotically flat spacetimes: 

Given "reasonable" asymptotically flat initial data of Einstein-matter field equations and 
assuming "reasonable" energy conditions, there exists a Cauchy evolution that is global in 
the sense, that a solution does exist outside black holes for arbitrarily large values of the time 
of an asymptotic observer . 

It is possible to formulate a " dynamic" version of the cosmic censorship in cosmological 
models ||, where the geometric approach runs into a trouble because there exist serious 
conceptual problems in defining a notion of a black hole. Most of the existing literature is 
concerned with the validity of that version of the cosmic censorship in various cosmological 
models |J. 

An analysis of spherically symmetric systems suggests that the two approaches to the 
cosmic censorship can converge. Trapped surfaces, that are inherent to the study of the 
cosmic censorship conjecture by Krolak , have been investigated in spherically symmetric 
geometries HlQfl , |TT| . Their appearance is always connected with a large concentration of 



matter. If a sphere S that is centered around a point of symmetry is trapped, then it 
contains amount of matter M(V) = Jv(s) P^-V °f t ne order L(S) [fH]] , where V(S) is the 
volume inside S and is the geodesic radius of S. And conversely, a large amount 

of matter, M(V) > L(V), leads to the formation of an apparent horizon. The quantity 
fv(S) P^ can De bounded from above and from below by a suitable Sobolev norm (usually 
Hi) of a matter field. Sobolev norms are, in turn, natural objects in the dynamic approach; 
an evolving system does exist so long as its Sobolev norm is finite. Singularities can be 
defined as those regions of space-time that give an infinite contribution to a Sobolev norm. 
Take an initial configuration of a compact support of a hyperbolic system of matter fields 
coupled to Einstein equations. Hyperbolicity (I always use the Petrovsky meaning of that 
term; see Section 3) means that the support of matter remains finite during a finite evolution. 
Assume that a Cauchy evolution can be analyzed in the Sobolev class Hi. If an evolution 
breaks down then there exists a sphere S such that a contribution to the Sobolev norm Hi 
coming from the interior of S becomes big in comparison to L(S) - hence there must exist 
a trapped surface and an event horizon enclosing S. That means that, in the spherically 
symmetric case at least (and under the assumption that there exists an evolution in Hi - 
that fact is not obvious), the notion of strong singularity is related to the notion of dynamic 
singularity. 

The "dynamic" version of the weak cosmic censorship hypothesis leads to an external 
Cauchy problem with a free inner boundary being a null cone, as is explained below. Take 
a four-dimensional (asymptotically flat) Lorentzian manifold M and define its space-like 
foliation by hypersurfaces S t , where t is the asymptotic time, t < oo. A black hole H can be 
defined as the largest piece of M that still can be enclosed by a null cone SH such that the 
area of the intersection 5H H E 4 remains uniformly bounded. (That definition is motivated 
by the fact that an area of a spherical black hole is bounded from above by 16irm 2 , where m 
is the asymptotic mass |I2| .) The weak cosmic censorship hypothesis reads in this context 
as follows: for "reasonable" initial data and "reasonable" matter, there should exist a global 
Cauchy evolution of a region exterior to H. The inner boundary (event horizon) 6H of that 
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region is a null cone so that we arrive naturally to an external Cauchy problem with an 
event horizon SIS db free null- like inner boundary. 

In this paper I will study a more general case of a Cauchy evolution with a free null- 
like boundary that does not necessarily coincides with an event horizon. A selfgravitating 
nonlinear scalar field is used as a matter model. The main results of the paper are contained 
in Theorems 1 and 2 of Section 4. Their formulation would require a number of preliminary 
definitions, so instead let me just discuss the most important points. 

The fact that I find the most interesting is that even the local version of the Cauchy 
problem requires the positivity of the potential (selfinteraction) term W(4>) > 0. Using a 
technical jargon, negative selfinteraction W((f>) means that scalar field equations (see |28|, 
p9|) would loose their strict hyperbolicity. (^) implies that the loss of hyperbolicity is 
particularly plausible inside a region enclosed by the Schwarzschild radius 2m. 

The smoothness of an evolving solution is locally preserved, that is its differentiability 
properties are kept intact for a small enough interval of time. The global existence is also 
proven, through standard apriori estimates, in the differentiability class (peH^, (k > 3). 
Having the global existence outside a null cone SH, one is given a geometry of SH and a 
scalar field on SH; they in turn can be regarded as initial values of the characteristic Cauchy 
problem. The present work is supplementary to the investigation of Christodoulou ||13|| , who 
investigated the existence of a solution of the characteristic Cauchy problem for (massless) 
scalar fields. 

Although I deal only with scalar fields, my impression is that the formalism presented 
below is capable to include other forms of nontachyonic matter, that is assuming that matter 
fields are described by hyperbolic equations and their energy-momentum tensor satisfies some 
energy conditions. Spherically symmetric Yang - Mills and SU(2) nonlinear sigma model 
do in fact satisfy the needed conditions. Einstein - Vlasov system does not belong to that 
category that can be analyzed with techniques presented here, but let me mention that a 
recent work of Rein, Rendall and Schaeffer shows that there exists a global solution for 



Einstein - Vlasov system ||14|| , if matter does not form a singularity at the symmetry center. 



The content of the rest of the paper is following. Section 2 presents the Einstein equations 
coupled to a spherically symmetric field equations in (1+3) splitting. They can be reduced 
to a system of integro-differential equations with matter fields being dynamical variables. 
The reader may find it useful that one can explicitly express components of the spherically 
symmetric metric as some functionals of matter-related terms - see formulae (|nj - |TB|) - 
in any splitting of the space-time, that is for any value of the trace trK of the extrinsic 
curvature. 

Section 3 describes the Einstein - scalar field equations in the so-called polar gauge. That 
way of foliating a space-time is particularly convenient when discussing the external Cauchy 
problem with a free inner boundary. In most gauges one is confronted with the need to 
impose an additional boundary condition at the inner boundary; in contrast with that the 
polar gauge does not require any extra conditions. 

Section 4 contains the main result, Theorems 1 and 2. The proof of the local Cauchy 
solvability bases on the standard compactness method. The main part of the local part of 
the proof of Theorem 1 is relegated to Section 5 while the global part is based on important 
L p and Hi global estimates that are proven in Section 6. In Section 6 a change of dependent 
variables leads to a set of first order equations that are "almost" linear; that allows one to 
obtain a set of apriori global estimates. Section 7 discusses a generalization of the proof of 
Theorem 1 that finally leads to the proof of Theorem 2. 



2. Equations. 

The most general metric of a spherically symmetric spacetime is given by the expression 

ds 2 = -N 2 {r, t)dt 2 + a(r, t)dr 2 + R 2 {r, t)dQ 2 , (1) 

where t is a time coordinate, r is a radial coordinate, R is the areal radius and dfl 2 = 
d6 2 + sin 2 9d(j) 2 is the standard line element on the unit sphere with with the angle variables 
< <p < 2tt and < 9 < it. I assume that at spatial infinity N — 1 (hence t coincides with 
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the proper time of an observer who is far away from the center of symmetry) and a — 1. 

Initial data of Einstein-scalar fields equations on an exterior part E° u< of a Cauchy slice 
Et consist of (g^-, Ky, p, ji) where gij is the intrinsic metric, Kij = ^jjff- is the extrinsic 
curvature, p = — T® is the matter energy density and ji = NT® is the matter current 
density. Tpv is the energy - momentum tensor. Initial data satisfy the constraints 

(3) i? - KijK ij + {trKf = lQirp (2) 

ViK ij - VHrK = -8nj j (3) 

where is the scalar curvature of the intrinsic metric. Above (and elsewhere in the 
paper) I assume the Einstein summation rule with the exception of indices r, 6 and whose 
repetition is supposed never to mean summation. Latin indices change from 1 to 3 while 
Greek indices range from to 4. 

In the spherically symmetric case it is convenient to formulate the whole set of Einstein 
equations in terms of the extrinsic curvature (which describes the embedding of the 
three-dimensional hypersurface E into a spacetime) and the mean curvature p (which is 
a trace of the two-dimensional extrinsic curvature that describes the embedding of a two- 
dimensional sphere centered around the symmetry center into E t ). The interior metric of 
E 4 is 

ds 2 {3) = adr 2 + R 2 dtt 2 . (4) 

The mean curvature of a centered two-sphere as embedded in an initial three dimensional 
hypersuface is 

2d r R 

If the trace trK of is fixed then there is only one independent component of the 
extrinsic curvature, say the radial-radial component K = and the remaining two com- 
ponents are each equal to *' h " A 



2 ' 
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trK-K = 2K 6 e = 2K$. (6) 
In terms of K and p the constraints can be written as 

J^J. = _ 8vrj R p - —{ K f + T (trK) 2 + -KtrK - -f- + - (7) 

and 

d r (R 3 (K -trK)) , j r r ,„„ 

— = -8ttR 3 ^= - ptrKR? (8) 

In the spherically symmetric case the full set of the Einstein equations consists of the two 
preceding ones, the evolution equation 

d (K - trK) = ^-(K) 2 + ^(trK) 2 - 2NKtrK - ^d r ^- + 8nN{T r r + p) (9) 
2 2 y 1 a pR 

and the lapse equation 

V^N = N(^{K) 2 + {trK) 2 j2 - KtrK + 4n{p + 7?)) + d trK. (10) 

The Einstein equations and Bianchi identities yield the energy-momentum conservation 

equations, which in the case of spherical symmetry reduce to 

i i N r) N 

d ^= + N(K + trK)^= + ^d r T r r + ^(T r r + p) + Np(T r r - T*) = (11) 

- d p - *Lj r - ^d r ^= - ™!*j r - NK(T r r - Tp) - NtrK(p + T+) = 0. (12) 

Using the above equations, one can express metric coefficients in terms of T^ v and trK. 
From the momentum constraint one finds 

„ C + 8n J™ 3r r 3 dr J™trKr 2 dr 
RK(R) — RtrK(R) = j^- + 2^-^ . (13) 

The parameter C is constant on a particular Cauchy slice and it must be set to on slices 
including the symmetry center. C is arbitrary, however, on slicings that do not include the 
world line R = 0. 

The integration of the hamiltonian constraints leads, after some algebra, to the expression 
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? R = T - f + ^ + T (frf " A '» *• (14) 

where m is the asymptotic (ADM) mass. The function m(R) can be interpreted as a local 
energy energy density (it is easy to notice that m(0) = m) and it is given by the equation 

m(R)=A7T drr 2 (p + rj r (trK — K)). (15) 

J R 

The lapse N can be determined from 

w= g (l+4 £» *cg^^ )w (16) 

where 

0( r ) = jr^- T r-p)- 2KtrK + 2 ( trK ) 2 )^r s d \ ( 17 ) 

The lapse N satisfies (|K]), which can be shown by using the conservation equation (fLip . 
The line element can be written directly in terms of p, K, trK, R and N: 

ds 2 = dt 2 (-N 2 + N2 ( trKR - KR ) 2 ) _ 2N trK ~ K dtdR + -4^dR 2 + R 2 dQ 2 . (18) 

(pR) 2 p [pRy 

Equations (0 - [18]) demonstrate that all metric functions can be expressed as certain func- 
tional of matter-related terms. That is simply a manifestation of the well known fact that 
spherically symmetric gravitation does not carry degrees of freedom independent of matter. 
The scalar field equation can be cast into one of two equivalent forms 

—6' 6' d N 

DD<p = —D(f)- -D<p- j^D<f)-trKD<f) + W', (19) 

DD(j) = ^-Dcf) - °-D<p - jjj=D<j> + trKDcj) + W. (20) 
Here the differential operators D,D_ are defined as follows 

£>=tU + -U, £=~ft + -^ (21) 



6' = p + K — trK, 9 = p — K + trK are the optical scalars and W = d^W{4>). W(4>) is the 
scalar field selfinteraction potential. 

8 



The nonzero components of the energy-momentum tensor of the scalar field are 
T ° = ^( {D<pf + GD0) 2 + r r = \{{D<pf + - AW{<j>)) 

% = w im)2 ~ ( - 0)2) ' T * = T * = ~Y wm ~ 2wm - 



(22) 



3. Definition of the external Cauchy problem. 

We will deal with an external Cauchy problem with a free inner boundary. Define £q u * 
as an open end extending outside a sphere of a coordinate radius Rq. 

Initial data entirely determine the geometry of the initial slice and, given the conservation 
equation of the energy momentum tensor, the lapse function N; see equations (|TB|- ^). That 



specifies also null cones attached at E . Take now an outgoing null cone H originating at 
a radius Ro. The open end £q m * shall give rise to a foliation defined by open Cauchy ends 
E° ut ; notice, however, that the foliation is nonunique in general. That is because on all 
future (t > 0) open ends E°"* the extrinsic curvature K (see ([13])) depends on a parametr C 
that is arbitrary. One would have to impose a new condition on the free inner boundary (in 
addition to the fixing of the trace of extrinsic curvature trK of ££"*) that should guarantee 
that the lapse function iV is strictly positive in H° ut . This is in order to guarantee that the 
proper time r = / Ndt runs forward in the outer region. This demand is quite restrictive, 
in general. ^From formula ([16]) follows that the lapse can vanish on a centered sphere S if 
the mean curvature p(S) = 0, but it might vanish also elsewhere, if the second factor of ([IB]) 
equals 0. The second posibility is more difficult to deal with, since the factor in question 
depends on the acceleration of matter (see ( |TTD and (0)). Fortunately, the arising difficulty 



can be avoided by imposing the polar gauge condition trK = K\ that condition removes 
the arbitrariness that is present in other gauges so that there is no free data at the inner 
boundary. Indeed, in such a case equations (]14| -|I8|) become 



ri „ r 2^ 2m(R) 

pfl=2 V 1 -if + ^ i ' (23) 



("OO 

m(i?) = 4vr / drr 2 p, (24) 



N = ^(3(R), (25) 



/3(r) = e 167r r(-^-")^ ds . (26) 



and 



ds 2 = -dt 2 N 2 + —^-—dR 2 + R 2 dQ 2 . (27) 
ypRf 

The parameter m appearing above is a free parameter that will be assumed to satisfy the 
inequality m > m(_R ); the equality might take place only if the energy density p vanishes 
at the inner boundary (exact conditions are given in Theorems 1 and 2). ( ]2"3"| ) - ( [ZED imply 
that N is positive provided that p > 0. The momentum constraint @ can be used in order 
to replace trK by = ^ L j r / \/a and the equation @ allows one to eliminate the gradient 9 r iV. 
With all that equations (19) and (|20|) lead to 



(ft + iMfcjy = - T) - ^(U + K/2) - ^ + W W , (28) 

(fib - ^ffiW = ^fi/0 + 71 + f (F + U/2) + ™ - W(«; (29) 



above we set 



y = £0 

U = D(f) 

jr 

'a 



rp rpr 

J- r • 



(30) 
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We will refer to fl28|) and (^) as to "reduced Einstein-scalar field equations". The whole 
dynamics of the selfgravitating scalar field is given by the scalar field equations (|28| ) and (p5| ) . 
The scalar field <p can be written as an integral quantity 4>(R, t) = 4>(R, 0) + J * drN(U — V)/2 
so that functions U and V can be used as dynamical variables. The initial data of the whole 
system of Einstein - scalar field equations consist of U and V supplemented by a value of 
at a single point (say 0(t = 0, oo) = 0). 

Equations (|8|) and (|9|) are strictly hyperbolic in the sense of Petrovsky provided that 
NpR is strictly positive. One can show that NpR can vanish only if pR = 0. On the other 
hand, d r /3 (and hence some of coefficients of the reduced equations) can become singular at 
those points, where the mean curvature p vanishes. Thus in our case strict hyperbolicity is 
a prerequisite for the existence of a local causal evolution. That means, from the inspection 



of (|23] and (p5]), that close to the sphere with R = 2m the contribution m(R) has to be 
positive, that is guaranteed if p > 0; the latter condition is satisfied if the selfinteraction 
W(4>) of the scalar field is nonnegative. Thus, even the local version of the Cauchy problem 
might require the positivity of the potential term W(<f>) > 0. 

As pointed above, the system of equations fl2%D and (^) is strictly hyperbolic in a region 
that does not contain minimal surfaces. From the inspection of formulae for p and N it 
is obvious that both quantities are strongly positive outside a sphere of the Schwarzschild 
radius R = 2m, so that outside that region the information propagates causally and time runs 
forward. That would mean that the sort of the external Cauchy problem that is described 
in the beginning of this Section is legitimate, at least for initial open ends Sq"' originating 
out of the sphere R = 2m. We will show that to be true even in a more general situation. 
We shall stress, however, that the coordinates that we use do not allow for the investigation 
of regions inside apparent horizons because in the polar gauge apparent horizons coincide 
with minimal surfaces, where our system of coordinates breaks down. 
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4. The Cauchy solution. 



The reduced equations f29| ) are integro-differential and hyperbolic; notice that they 
are nonlocal. I did not find any mathematical result that can yield directly the existence 
of a local in time solution for equations of that type. For that reason we have to prove the 
existence of a local evolution starting from first principles [R] . 

There are several methods to prove the existence of a local Cauchy solution. We will use 
an approach that bases on results of Petrovsky || and on properties of Sobolev spaces. The 
main result is formulated in Theorem 1. Its proof consists of following main points. Firstly, 
a sequence of functions will be generated iteratively for data given on an extended initial 
hypersurface. Secondly, Lemma 3, proven in the next Section, shows that the sequence is 
uniformly bounded for a short period of time outside a "rigid" cone. Then, in Lemma 4, 
standard compactness theorems of functional analysis ensure the existence of a convergent 
subsequence, that is the sought local in time solution of the reduced equations. Finally, 
Lemma 7 shows the existence of a global Cauchy evolution. 

Let us recall that Sobolev spaces H k (V) can be defined as a completion of C fc -functions 
in the norm \\f\\ Hk (v) = Jy dVEf =0 (D l f) 2 where D % f = Z klM ... kn . kl+k2+ ... kn=i d kl ...d k J and 
n is the dimension of a riemannian manifold V. In the case of spherical symmetry the norm 
reads ||/|| = J c?rr 2 Sf =0 (9*/) 2 . The main result is the following one. 

Theorem 1. Let the initial data of the reduced equations (|2^ ) and ( f2~9|) be U,VeH k 
and <peH k+1 R L no , k=2, 3, ... . Assume that W{x) > for x ^ and W(0) = and 

\d l x W{x)\ < A{l,m)\x\ n ' 

for some < n\ < oo and a constant A depending only on n\ and I = 0, k. Let Sq u * be an 
open end, Rq = inf [|x| : xeSg"*] > 2m + rj, rj > and m > m(Ro) be an asymptotic mass of 
the configuration (with strict inequality m > m(R ) if d l R f ^ 0, / = U, V, at least for one 
value i — 0, ...k — 1). 
Then 

i) there exists a local Cauchy evolution of Sq"*, i. e. a foliation E° u * for some r]'( 

12 



< t < rf), with U, VeH k and <peH k+l . In addition, U, VeC k _- l ({j Q < t<rf S t ou '); 

ii) solutions are unique; 

iii) solutions are global in H k and the null boundary 5H of Uo<t<r/ ^° ut escapes to spatial 
infinity, i. e., the area of the inner boundary goes to infinity as t —>■ oo. 

Actually, one can prove a sligthly stronger version: 

Theorem 2. Let the initial data of the reduced equations (^8| ) be U, VeH k and 4>eHk+in 
L no , k=2, 3,.... Assume that W(x) > for x ^ and W(0) = and 

\d l x W(x)\ < A(l,m)\x\ ni 

for some < rii < oo and a constant A depending only on n\ and I = 0, k. Let Sq u * be an 
open end, Rq = inf [|x| : xeSg"'] and m > m(Ro) be an asymptotic mass of the configuration 
(with strict inequality m > m(R ) if d l R f ^ 0, / = U, V, at least for one value i = 0, ...k — 1). 
Assume that minimal surfaces are absent in Eg"'. Then 

i) there exists a local Cauchy evolution of Sg u *, i. e. a foliation S° u * for some r]'( 
< t < V), with U, VtH k and <ptH k+l . In addition, U, VeCt^iUoKKr,' ^7% 

ii) solutions are unique; 

iii) solutions are global in H k and the null boundary 5H of Uo<t<r/ ^° ut either 

a) escapes to spatial infinity, i. e., the area of the inner boundary Y>° ut goes to infinity 
as t — > oo, 

or 

b) stabilizes at the radius R = 2rriB, where tub is the Bondi mass of a black hole. 
Comments. The condition that m may be equal to m(R ) only if a sufficient number of 

derivatives of U and V vanishes, is always satisfied if initial data on Sq"' are obtained by re- 
striction from initial data on a whole slice (including the origin) S . Conditions U, VeH k (T,Q) 
mean (through Sobolev embeddings theorems) that k — 1 derivatives of U and V are contin- 
uous. If m = m(Ro) then U(R) and V(R) must identically vanish for R < Rq, and therefore 
also their derivatives up to the order k — 1. The symbol C*_ means a class of functions which 
are almost C k , that is they are of Holder class C k ~ 1+K for any < k < 1. Let us point 
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out that the smoothness of U and V can be improved to C^ 1 ^ 2 , at least for initial data of 
compact support, following an argument of ||15|| . 

In what follows we will prove Theorem 1 in four steps and five lemmae; the proof of two 
of them is placed in separate sections. 

Theorem 2 can be obtained by improving one of the elements in the reasoning. 

Proof. Step I of Theorem 1. 

Define a extension Eq u * = [xeS : 2m+i] — 5 < |x| < oo] of Sq"* by assuming following 
initial data in the interval [2m + rj — 5, 2m + rf\ 

d k J(R) = 0, d^f(R) = d^fiRo), d k T ~ l f{R) = #~7(i2o) - / drd^ +1 f(r), 

JR 

where f = U,V and i = 2, ...k. The initial values of the scalar field in the extension can 
be obtained from integrating U and V, 4>{R) = 4>(Rq) — Jr° drj^(U + V). One can check 
(see Appendix) that norms of U and V over the extended manifold Sq u * differ from 
norms of U and V over Sq u * by a term bounded above by 8CH; t (T,Q it ), where C is a constant. 
Similarly, the mass function m(R) < m(Ro) + <5C-£/i(£q u *) for Re[2m +rj — 5, 2m + 77]. Thus 
for S sufficiently small the extended manifold can have the same asymptotic mass m and lie 
outside the Schwarzschild radius 2m. 
Let us define following quantities 



P0 (R,t) : 


= P (R,t = 


0) 


N (R,t) 


= N(R,t 


= 0) 


U (R,t) 


= U(R,t = 


= 0), 


V (R,t) 


= V(R,t = 


= 0), 




= <f>(R,t = 


= 0) 



as being identically equal to the initial data. Let po, To and jo be the initial values of the 
energy density p, T and j, respectively. Define also, for n > 1, functions 

f 2m-8n fx drr 2 p n (r,t) 
a n {K, t) — y 1 — 
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-4tt f°° dvr {Tn+ , Pn) 

Pn = \((U n f + (V n f + 4W(<J) n )), T n = \((U n f + (V n ) 2 - AW(<p n )) 

Jn = \{{V n f - (U n ) 2 ) } 

<f> n (t, R) = <f>(R, 0) + J* dra n _i(r, R)P n -i(r, R){u n -i(r, R) - K-i(r, R))/2. (32) 

All formulae above are expressed in terms of U n and V n which (for n > 0) are defined as 
solutions of following evolution equations 

(d + u 2 n -il3n-xdR)V n = 

47rP n ^RV n (j n - T n ) - Q "-^ n ~ 1 (U n + V n /2) - + a n ^ n ^W'(<f> n ), (33) 



(d - al^Pn^dnjUn = 

47r^xRU n (j n + T n ) + a ^ n -\ V n + U n /2) + %^ - a^x^W (cj> n ) , (34) 

with initial conditions U n (R,0) = U (R) and V n (R,0) = V (R). 

Define tt t = [(t',x) : < t' < t, |x| > 2m + r] - 5 + 2t'] and £f* = [xeVt t : x° = t]. 
Define 

out \ / (I afc r r 1 2 i I ofc t /" 1 2 



£UEr*) = / i^xr + KVnV). (35) 

The essential part of step I is following. 

Lemma 3. Let the initial data of the reduced equations (^ ) and fl29|) be of compact 
support and U , VoeH k (Y, Q Ut ), k > 2), for k=2, 3,... Assume that 

W(x)>0, \d l x W(x)\ < A(l,m)\x\ ni \ 

for some < n\ < oo and a constant A depending only on n\ and I. Then, for a small 
enough time t: 



i) for each n > there exists a solution (£7 n , V n ) of equations ( p3| , |34l ) in Q t ; 

ii) Sobolev norms H k (Y,° ut ) of U n and 14 are uniformly bounded for t' < t, 
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EmVT) < 7 rim*, (36) 

where I = 0,1, ...,k, E = EiniT^ 1 ) is the initial value of the norm E\ n and CVs are some 
constants; 

iii) coefficients a n ,(3 n ,(j) n of ( |33| , |33|) are at least C k as functions of (t,R) with a norm 
l/lc* = sn P^k 1 +k 2 =o\^o 1 ^R f\ (where / = ot n , /3 n ) uniformly bounded by (up to a constant 
coefficient depending only on initial data) 

K(sr)f < — — 1 , 2fc/(fc+2) (37) 



and, in particular, a 2 /3 ra < 2; 

iv) L 2 (S° U< ) norms of d l Q d^f (f = Ui, V\ and I > 0, / + m < k) are uniformly bounded; 

v) \d l R (p n (t, R)\ < C^ + C^t, where C3 and C4 are constants depending only on initial data 
and < / < k. 

We postpone the proof of Lemma 3 to the next Section. Assuming that it holds true one 
shows the following 

Lemma 4. Assume conditions of Theorem 1 and, in addition, that initial data are of 
compact support and k > 2. There exists a solution (U, V)eHk(D, t ) and (U, V)eC k ~ 1 (Q t ) of 
the reduced equations (p8| , ^) in fl t . 

Proof of Lemma 4. Notice that there exists a compact set Q t that contains supports of 
all functions (U n , V n ). ^From iii) and iv) of Lemma 3, U n , V n constitute a bounded subset of 
Hk(Qt)- Therefore there exists a subsequence U nk ,V nk that is weakly convergent in Hk(flt) 
and strongly convergent in iffc_i(f2 f ) to a limit (U, V). (U,V) is H:older C M continuous, 
hence (U, V) is a limit also in the norm. (U, V) solves, by the construction of the 



iteration, the reduced equations (|28], |29|) in Vt t . The solution belongs to H k {Vt t ) and, since 
our problem is essentially two-dimensional, it is C k ~ 1+K (Qt), for any < k < 1, by the 
Sobolev embedding theorem. That ends the proof of Lemma 4. 
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Lemma 5. Under conditions of Theorem 1 and for k > 2, a solution (U, V) is unique in 

fit. 

The proof goes in a completely standard way and we omit it. 



Step II of the proof of Theorem 1. 

We can remove the compactness condition of the support of initial data. 
Lemma 6. Under conditions of Theorem 1 a solution (U, V) with (noncompact) initial 
data that are of class Hk, k= 2, 3,..., exists in Q t . 
Sketch of the proof of Lemma 6. 

We shall discuss the case with k = 2 only. Take a C°° sequence of functions of compact 
support U 0n , V 0n that approximates (U , V ) in H 2 (T>Q Ut ). Each pair of the sequence U 0n , V 0n 
gives rise to a solution U tn , V tn of the reduced equations (^, |29[) in Vt t - From the reduced 
equations one derives, after a lengthy but simple algebra, 



Thus the sequence of solutions is uniformly bounded in H 2 (flt), for sufficiently small but 
nonzero t > 0. Therefore there exists a subsequence (U tnh , V tnk ) that is weakly convergent to 
(U, V) in H 2 (tt t ) and strongly convergent in H x (nf), where Qf = [(R, t)eQ t : R < M]. The 
last fact and the Holder continuity of (U, V) and (U tnk , V tnk ) imply pointwise convergence of 
the subsequence in Qf 1 . A careful analysis of the reduced equation shows that for M » 1 
the reduced equations are nearly linear. A perturbation analysis shall prove the existence of 
a solution in Q t \Qf ! that matches to a solution (U, V). That concludes the proof of Lemma 



Step III of the proof of Theorem 1. 

The inner boundary 5Q t of &>t consists of pairs (t, 2m + r] — 5 + 2t). The outgoing null 





3 



(38) 



6. 



rays in Vt t are determined from the differential equation 



/From (H - H) 



wc 
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conclude that ^ < 1 at 5Q t . From the smoothness condition we infer that l^L^ 1 / 2 [ s a t 



least a bounded function along a null ray. Thus 5Q t is a space-like three-surface. Therefore 
there exists a null outgoing geodesic joining each point (t, R) (for small enough t) of Sfl t 
with a point P on the initial open end Eg"'. Let us choose t such that the point P lies on 
the boundary of Eg"'; that * s possible since by definition Eg"' C Eg"'. Then there exists a 
time development E°"' C E°"', defined by the solution (U, V) of the reduced Einstein-scalar 
field equations, of the initial open end Eg"', i. e. Eg"' is future null geodesically complete. 
That accomplishes the local proof of Theorem 1. 

Step IV of the proof of Theorem 1. 

Lemma 7. Define g = ^ and h = Under conditions of Theorem 1, 
Proposition 7.1. Banach norms L 2 j of g and h are equivalent to norms L 2 j of U and V, 
respectively; 

Proposition 7.2. Sobolev norms iJ; and integral norms L 2 j of g and /i are uniformly 
bounded by Ce ct for any integer j, I < k and for some constant C; 

Proposition 7.3. The scalar field is uniformly bounded by a constant C depending only 
on an asymptotic (ADM) mass of the configuration and on the areal radius R of the initial 



The proof of Lemma 7 is highly technical and we postpone it to Section 6. Using the above 
Propositions we infer that, under conditions of Theorem 1, there exists a uniform bound on 
L 2 j norms of U and V. We show that norms of U and V remain bounded during a finite 
evolution. Indeed, 



slice Eg"', 



\4>{R)\<C. 



(39) 




(40) 
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for some C. The first inequality is just a trivial application of 2ab < a 2 + b 2 and in the second 
one we use < pr < 2, expression fl23| ) (with p = \{U 2 + V 2 ) + W(4>)) and Proposition 
7.3. The last expression of (PH) is globally bounded by applying the Schwartz inequality and 
then using Propositions 7.1 and 7.2. Thus the norm | \d r U\ || 2 ( St ) is globally bounded. 

Similarly one shows the global boundedness of ||d r V||| 2 ( Et ) and, iteratively, of higher 
Sobolev norms of functions U and V . 

One can show that the area of the inner boundary of E t goes to infinity as t — > oo. 
Indeed, the inner boundary [j t E t is null and its areal radius Rt satisfies the differential 
equation ^ = ^ pR ^ (3. Initially we have Rq > 2m + i], with r\ > 0; thus from (|23D, (^) and 

~ 47r f drr 2 (p+T) 

the equation we infer that (pR) 2 (t) > 4(1 - and (3{R t ) > e 2m(1 -^ } Rt > 



; 2(1 ^Tv\ Thus is strictly positive nad R t goes to infinity. That completes the global 



part of Theorem 1. 



5. Proof of Lemma 3. 



The approximating equations ( p3|) and (134]) are nonlinear in U n ,V n , with polynomial 
nonlinearities in each order n, so that their coefficients are of class C°° as functions of 
unknown variables U n ,V n . A long but simple calculation yields following equations 

„2 

//>■ 2 



dr-((d k R U n ) 2 + (d k R V n ) 2 ) 



R 2 



R 



((d k R U n ) 2 (-2 - e^-i/Vi) + (4K) 2 (-2 + aLx^-x)) 



drr' 



+ 



r n - ^ ^K-x/^-i) + fc9 r (<_A-i)J ((<9X) - (d k v n y) + 



]7irp n . lJn (d k V n ) 2 + d k r {an-iPn-iV'){d k r V n - d k U n ) 



+ 



1^ drr 2 ^ =l [(dir n Q + (i -<yd;+V^A-i)(* + J)^ + 



(41) 



where 
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T n = 47r(3 n ^R(j n + T n ) + 



C^n—lPn—l fin 



+ 9.R ' U 



fc! 



2R ' 2i? ' Vi y i!(jfc-i)!" 
In particular one obtains 

d / r 2 dr Pn = -iyl - Ul)al_ x p n ^R 2 + 
Jr 4 

- y rfrr 2 /3 n _! [W(0 n ) (K - K-i -U n + Un-i) + l&Krj n (W(<j> n ) - W{<j> n 
Let us recall an obvious inequality 



(42) 



(43) 



1/2 



that holds true for any function of compact support that belongs to H\. 
One gets from (fill) , applying the Schwartz inequality and pointwise estimates 
inequality 



(44) 



the 



d 



dt 

fc-i 



sup |r n | + supp„,_i + sup \j n \ + J E^ 2 ||^(a„_i/5 ri _iV / )IU2 + 

£ ^K 2 (su P |#r„| + sup i^+Vn-A-i)!) + 

i=l 

< 2 < 2 (ll<tr„IU 2 + lie^atA-i)!!^ 



(45) 

Here and below constants may change from a line to line, but they do depend only on initial 
data and k. Using (fHj) and various Holder estimates in order to bound the sup — terms on 
the left hand side of fl4"5|) one derives the differential inequality 



{dt + 2d R )E kn {Y,™ t ) < C 2 E kn (E 2 ln + + E k k 



kn—l 



(46) 



The interesting fact is that for k > 1 the inequality ( f4"6"|) is linear in the term E kn . 

I will use the method of exact induction in order to prove Lemma 3. For n — 1 the 
coefficients of (^) and are time- independent. Sobolev embeddings theorem (or, strictly 
saying, (pH))) implies that Uq and Vq are of Holder class C fc_1 / 2 . Therefore the coefficients of 
(|32|) are of class (7 fc+1 / 2 for i? > Rq as functions of i? and (being time- independent) they are 
bounded and continuous in the class C k , k > 2 as functions of t, R. The eigenvalue «o/3o of 
the hyperbolic operator is bounded from above by 1 and it is strictly nonzero so that in the 
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case n = 1 the approximating equation is hyperbolic. A result of Petrovsky || ensures the 
existence of Ui, VieC k ~ l (Q t ), confirming part i) of lemma 3. 
16f) becomes now 



(d t + 2d R )E kl (Z™ t ) < C 2 E kl [E 2 n + E lfi + E*_ 10 (E k .. n + E k0 )). (47) 

where E k0 = E^iY*™ 1 ) is the initial value of the norm E kn and Ci,C% are some constants 
depending only on initial data. The C k ~ x smoothness of the solution U\, V\ guarantees that 
En is a continuous function of time t for / < k. Now, let Ea(t) > Eq for t' > t > 0. Then 
([I?]) implies, after elementary calculations, 

^ * (E m '-> -C4" M ' ' <48) 
where we used the initial value condition En = Eiq. The above trick does not work for E k \ 
since the latter need not be continuous. In order to show a corresponding estimation for E k i 
one has to insert estimations fl4"8|) into fl4"TD and integrate the latter. The linearity of (|7 



in 



Eki allows one to get an estimation as in (081). That proves the point ii) of Lemma 3 



In order to show iii) let us notice that d R a n or d R (3 n contain derivatives of U n and V n of 
the order k — 1 at most; that is functions ot\ and (3\ are of class C h for fixed t. That can be 
shown quite generally, in any order. A direct calculation shows that 

d a 2 n = —d / drr 2 p n (49) 
a j r 

(fi9|) and ( |43"1) show that the order of differentiability of <9oa n is equal to the order of dif- 
ferentiability of U n and V n , if coefficients of the preceding order, (3 n -\, a n _i are of class C k . 
Therefore d ai is C fc_1//2 in Ti° ut and C k ~ x in fi t because Z7i, V\tC k ~ x l 2 and they are C fc_1 in 
Qf One shows that 



/■°° r Zo 

d /3 n = -An(3 n / drr 2 -J^d a 2 n + 
Jr 1 rat 



2j n a 2 n _ 1 l3 n -id r — ? 

red 



1 



A,- 



ret 2 



- C/ n ) + 16vrr Jn (^(0 n ) - 
21 



71—1, 



+ 



-8Tr(3 n j n R 2 

o Ci n —\Pn- 
\/4 



obviously doP% is C fc ~ 1//2 in Thus both «i and /?i are at least C k in That accom- 

plishes the first part of the statement iii). C k ~ l - differentiability of Uq and Vq implies that 
also 4>i(t, R) = 4>o(R) + Jq dsaof3o U °~ v ° is (at least) C k . The estimate v) is trivially true for 
<pi. The second part of iii), the estimate (|37|) and that of afPi, can be obtained from C k 
estimates of a>i and Pi, expressing them in terms of initial data (through (f48|)). 

Differentiation of the approximating equations (|33| , [34|) with respect t and R, integration 
over Tj° ut and the use of ii) and iii) allows one to express those norms of functions Ui and 
Vi that are specified in point iv) in terms of initial data. That finishes the proof of the first 
step of the induction hypothesis. 

Now, let it be true for some n. Using the induction hypothesis and the same reasoning 
as above, one shows that the coefficients of the n + 1 equation ^) are at least C k , 



so that there exists a solution U n+ i and V n+ i, by a result of Pietrovski ( |5|). The second 
power of its Sobolev norm, E^ n+ i (where I < k), is either less than E^ n for some t close to 
t = or greater than E^ n . In the former case it satisfies the bound of (|48|) by the induction 
hypothesis while in the latter case we have, from (|46|) 

(d t + 2d R )E l , n+1 (f% a ) < {E lin+ i{t° t ut )) l+ \ (50) 

which again yields 



E hn+ i(Y>° t ut ) < — 1/(/+2) - (51) 

Eit 2 ~ Qt] 



The estimation of Ek n +i can now be obtained in the same token as that of Eki- Hence i) 
and ii) are proven. The remaining steps, iii) and iv) are shown identically as in the case 
n = 1. We will check the validity of v). Notice that 

0n+l = 00 (R) + / dsa n P n 



- - it,}- ib n 5 

o 2 



t 1 /2 

the right hand side of that equation is bounded from above by sup \4>q(R) \ + C j Q dsE ln , 



using fl44|) . The induction hypothesis ii) and direct integration immediately yield v). That 
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ends the proof of Lemma 3. Let us point out that the uniform bounds proven above show 
the existence of all approximating functions in a strictly positive time t; from (|51|) follows 
that t can be bounded from below by a number that depends only on initial data and does 
not depend on n. 



6. Proof of Lemma 7. 

Let R t = inf E cmi(r); we have 

r* NvR r f NvR 

R t = R + ds{d t + -^—d R )R =R + ds^— > R > 2m (52) 
Jo 2 Jo 2 

(these inequalities follow from the assumption made in Theorem 1 that Sq m * is placed outside 



a centered sphere of the Schwarzschild radius 2m). Using (14) and 52), we conclude that 

(53) 



2m 



Rp(t,R)> 2^1- — >0 



on all future slices T,° ut . Since (14) implies also Rp < 2, we infer that there exist two nonzero 
and finite numbers a, (3 such that 

V U 
otg < — < /3g, ah< — < f3h; 
pR pR 

that fact trivially implies the equivalence of L p norms of pairs (g, V) and (h, U), for any 
p, and hence also for even values of p. That proves Proposition 7.1. 
One can show that 



d (pR) = 87rNRj 



(54) 



The reduced equations (p8|) and (p9|) lead now to following system of equation for functions 



9 = aIld h = §R 



fr . NpR^ . ic Ng , . 2N Nph N TX/I ,,. 
(do + -^d R )g = 167r-^(0) -— g -^ + _^), 



NpR 
( d o i^d R )h 



Nh , s 2N , Npg N . 
-^-^(0) + ^ + ^-^(0). 



(55) 



(56) 
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The interesting fact is that these equations are almost linear in unknowns; although coef- 
ficients are some functionals of g and h, they are pointwise bounded by constants that are 
independent of both functions. 

Embeddings theorems imply that initial data (U, V) that are of Sobolev class Hk (k > 1) 
must belong to L p {Tj° ut ) for any p > 2; that implies, taking into account the definition of g, h 
that initial data of the the latter are also of class L p . Let n be an even number. Multiply 
(IH) and flSlf) by g n ~ x and /i n_1 , respectively and integrate over E° u *. Integrating by parts 
and using ( |61~D (see below) one arrives at 

?3 



(ao+ ^ R) r^v+ 9 ") = -«"+ 

2 Jr n n 



drr 2 



R 



pr v 4nr nr n 



9 Af(ft"- 2 - 9 "- 2 ) + ^^(j"- 1 -ft"- 1 )]. (57) 

The assumptions of Lemma 7 (a;) > for i^O and Vt^(O) = imply that for initial data 
of compact support the scalar field <fi vanishes at spatial infinity; then one can use fl4*4] ) and 
replace d r (j) by r^(£/ + V) (see ([3(]), which leads finally to the estimation 

<p(R,t) < c(\\u\\ L2 ^r) + \\v\\ L2{ ^ t) ) < c(\\u\\ L2 ^r) + \\ v W^rY> ( 58 ) 

the last inequality is valid only if the selfinteraction term W(4>) is nonnegative and it follows 
essentially from the conservation of asymptotic mass m. The constant C can be found 



explicitly and the right hand side of (|58|) can be found to be bounded by , = yl^- With 
this we prove Proposition 7.3. 

With the inequality (|5^) and having hitherto known estimate pR < 2, we get from 



(d + ^d R ) r drr 2 \h n + g n )<C f°° drr 2 -(h n + g n ) } (59) 
2 Jr n Jr n 

where C is a constant depending only on n and initial data. Therefore L n norms of g and 
h are bounded in a global time evolution. That proves the L p part of Proposition 7.2. 
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Differentiating equations (|55| ) and (|56[) with respect r, multiplying them by d r g and <9 r /z, 
respectively, and integrating over H° ut one arrives at 



d 



outj R drr 2 -((d r g) 2 + (d r hf 



dt Jr 

, >3 (<9 r /*) 2 + r drr 2 \({d r gf - (d r h) 2 ) (=±d r (Npr) + 16vr^ -— 2 + 
2 Jr ' v 1 v vr z 2 



NpR 3 



d r (-W'{<t>))(d r g - d r h) + (gd r g - hd r h)d r {16ir— - ^) + \{gd r h - ^)<9 r (iV>)j (60) 

Notice that coefficients ^d r (Npr), Wn ^ w _)_ ^ j n front of (d r h) 2 and (d r g) 2 are pointwise 
bounded on E° ui . That follows from obvious estimates < N,pR/2 <1,R> 2m, from the 
equation 

-^d r {Npr- 1 ) = Z ^r + - 47riVp- 1 r(T - p) = 

-_^V + 3A^r + g7rNrW ^y p ( 61 ) 

and from Proposition 7.3. Using the Schwartz inequality, Proposition 7.3 and the already 
proven part of Proposition 7.2 one obtains from (|60| ) the inequality 



✓7 /"OO 1 , . roc 1 , . 

dt |out /* ^2 ( {8r9)2 + (9r/l)2 ) " C /e dTr \ ( {dr9)2 + {drk)2 ) (62) 



with a constant C depending smoothly on various norms L^j and on the estimation of the 
scalar field. From (^) we infer that the integral drr 2 ^ ((d r g) 2 + (d r h) 2 ^j remains bounded 
during a finite evolution. Therefore Hi norms of U and V remain bounded. 

Now, let there exists a global estimate for the Sobolev norm H^-i of g and h. In one 
dimension (to which our problem essentially reduces) the boundedness of Hk implies the 
existence of L p estimates of d l r f, f — V, U, for I < k and for any integer p > 2. Differentiating 
reduced equations fl55| ) and fl5"6|) k times with respect r, integrating the resulting equations, 
using various Holder inequalities and the induction hypothesis, one obtains an inequality 
analogous to (^) for the norm. That leads to the conclusion that norm is also 
globally bounded. That ends the proof of Proposition 7.2. 
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7. Proof of Theorem 2. 
Now I will discuss shortly the proof of Theorem 2. It proceeds in a way analogous to 
the proof of Theorem 1. In particular, as in the former case, one defines the approximation 
procedure of Step 1 of Theorem 1. The only difference between the two cases is due to the 
fact that previously we used an open initial end Eg"* with |x| > 2m while now it penetrates 
a region inside the Schwarzschild radius |x| = 2m. The difference is important, since the 
former choice guarantees strict hyperbolicity of all approximating equations (|33| , |34|) - all 



coeeficients a^/3 n are uniformly bounded away from - while now we have to show that. 

Let Rq < 2m be the areal radius of the inner boundary of Eg"*. Define Ri = Rq — r > 
and Eg"* = [xeE : R 1 < |x| < oo] (an open end containing Eg"*). Let S° ut = [xeE : R < 
r < |x| < oo] be an open end contained in Eg"*. Define a four-manifold Q t — [(f, x) : < 
t' < t, |x| > Ri + 2t'] foliated by hypersurfaces E°"* = [xeQ t ■ x° = t] and another four- 
manifold tt Srt = [(i',x) : < t' < t, |x| > r - 2t'\ with leaves S°f = [xtVL rt : x° = t]. 
Obviously, S° r f = S° r ut . 

As in Section 4 we define energy norms of solutions of the approximating equations 
33) for each family of foliations, 



Ekni&rt 1 ) — (i^r^nlx 2 (5out) + l^r ^nli 2 (S°««)) • (64) 

All iterative solutions satisfy the initial data, U n (t, R) = U (R), V n (t, R) = Vo(-R); norms 
of initial data are denoted simply as E^q and E^q respectively, for each of the two families 
of foliations. 

The main technical result of Section 4, Lemma 3, has to be replaced by 
Lemma 3'. Let the initial data of the approximating equations (|28] ) and ( p9|) be of 
compact support and let there exists a nonzero r such that 



inf 

Ro<r<oo 



Eg' - (2m - r + 2r) > C 6 t, (65) 
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where Cq is a constant defined below that depends on initial data. Let Uq, VoeHk(T,Q Ut ), 
k > 2), for k=2, 3,... Assume that W(x) > if x + and V{0) = and 



(\dlW(x)\<A{l,m)\x\ n ^ 



for some < n\ < oo and a constant A depending only on n\ and I. Then for t < t 
O) equations (|3~3"D and (0) are strictly hyperbolic in f^; 



i) for each n > there exists a solution (Z7 n , V" n ) of equations ( p3|, |J) in 

ii) Sobolev norms fffe(E^"*) of £/ n and are uniformly bounded, 

^n(Sr) < — rrmv , (66) 



where I = 0, 1, fc and Q's are some constants; 

iii) coefficients a n ,/3 n of ( |5B] , |3"1| ) are at least C k as functions of (t,R) with a norm 
\f\c k = sup Sj^ +fc qI^q 1 ^ 2 /! (where / = a n) (3 n ) uniformly bounded by (up to a constant 
coefficient depending only on initial data) 

( *^ r ^~-«r" - (67) 

iv) L 2 (S°"') norms of d l d^f (/ = £/„, 14, an d i > 0, Z + m < k) are uniformly bounded; 

v) |<9|j</> n (t, i?)| < C3 + C\t, where C3 and C4 are constants depending only on I and 
initial data and < I < k. 

Proof. The only new element in comparison with Lemma 3 is the strict hyperbolicity 
assertion, point O above. We have to show that the product a^/3 n is strongly positive for 
sufficiently small values of r. That will be achieved if we show that a n itself is positive, that 
is (from @) 

R + 2t-2m + 8n / drr 2 p n (r, t) > (68) 

JR+t 

for R> R\. That is obviously true for the n = 1 equation, since all functions with a suffix 
"0" (including a ) are time-independent and initially a is positive. Now let the n — th step 
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of the induction hypothesis be true. We shall show that the (n+ l)th equations are strictly 
hyperbolic. Manipulating with equations ( |33"D and fl3~3| ) one finds the equation 

/•CO /"GO 

{d -2d R )2j R drr 2 p n (r,t) = (d -2d R ) J r dr-((U n ) 2 + (V n ) 2 + 2W( ( j) n ) 



R 



^-{{U n ) 2 {2 - e^A-i) + (K) 2 (2 + «Li/3n-i)) 



+ 



<irr" 



(69) 



that describes the evolution of the n — th order external energy (modulo a coefficient 4tt) 
fx drr 2 p n (r,t) along the foliation formed by leaves Sfjf. /,From that we obtain (dropping 
out the positive boundary term and using pointwise estimates in order to bound the last 
integral of fl69|)) 

/•OO 

(d -2d R )2 drr 2 p n (r,t)> 

JR 

/ drr 2 87r/? n _irj n (p n - p n _i + T„_i - T n ) + V'a n _i^i(^i -U n + V n - K-i) > 

JR L J 

-C(sup(|j n | + |K| + | tf n | + |K-i| + \U n -i\)) jT° drr 2 p n (r,t), (70) 
where C is a constant. Using now ([0]) and (|35|) we observe that sup(|j n | + |V^| + \U n \ + 



|y n _i| + |C/ n _i| can be bounded from above by Ci£f n (£°f ) + C 2 £f n (£ 



1/2 



where C\ and 



C*2 are constants. Using now the induction assumption ( [SB] ) and using the obvious inequality 
y/x + x < 1 + 2x, we arrive at 

POO 

{d -2d R )2 drr 2 p n (r,t)> 

JR 



-c 4 (i 



(-^io 3 - Ci 



1/3 



drr 2 p n (r,t), 



(71) 



which yields 

/"CO poo 

/ drr 2 p n (r,i) > e~ Cbt / drr 2 p(r,t = 0) > 

poo 

/ drr 2 p(r, t = 0) - C 6 t. (72) 

J R 

All constants above depend only on initial values and the degree of nonlinearity of the 
potential term W{<p n ). The last constant enters the condition (p5T) stated in Lemma 3'. 
Using (|72l ) we can bound the left hand side of (^) as follows 
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oo 



R + 2t - 2m + 8tt / drr 2 p n (r } t) > 

JR+2t 



oo 

2 



R - 2m + 8vr / drrVr, t = 0) + t(2 - C 6 ). (73) 

JR+At 



Take now r of fl65D and notice that the condition (65) is satisfied by all values < t < r. 
Notice that E^ R+it = 8tt f™ +4t drr 2 p{r,t = 0); inserting the condition E^ R+4t - (2m - (R + 
4£) + 2t) > C^t > of fl65|) into ( ff3|) we obtain (r — £)(2 + C^), that is the last expression 
of ( [T3"D is positive for t < r. Thus we arrive at the desired result a^(i?, t) = i? + 2£ — 2m + 
87r fft +2t drr 2 p n (r,t) > 0. The approximating equations ([33|) and (|34| ) are strictly hyperbolic 
in f2 t . 

The rest of the proof of Lemma 3' is similar to that of Lemma 3. One can prove also 
results corresponding to those of lemmae 4-7, thus accomplishing the proof of a local part 
of Theorem 2 and the first half of its global part. I omit details. The proof of the remaining 
global statements of the point iii) of Theorem 2 is given elsewhere Jl5| . Below I will sketch 
the main idea. 

Assume that we have a complete initial Cauchy hypersurface with an apparent horizon, 
that is (using the polar gauge) with a minimal 2-surface on the initial slice. As a side 
remark let us point out that the area of an outermost apparent horizon cannot decrease (see, 
e. g., a proof outlined in |j"2"f); in fact it has to increase whenever matter crosses through 
the horizon, which moves acausally outwards. Asymptotically the spherically symmetric 
apparent horizont becomes null and its areal radius stabilizes at 2rriB < 2m, where m# 
is the Bondi mass. Thus all null geodesies starting outward from R > 2m# at t — are 
complete. The biggest set H that is still future null complete can be constructed as follows. 
Take a part S r of the initial hypersurface that does not include minimal surfaces. Then data 
on E r give rise to a local evolution, according to the local part of Theorem 2. The global 
evolution prolongs until the free inner boundary "freezes" near a minimal surface; the lapse 
collapses to values close to and the area of the inner boundary is practically constant. 
In such a case one can take a slighly smaller initial open end S r / C S r ; that evolves to a 
spacetime that freezes at a later time than the previous one. Continuing that procedure 
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ad infinitum one finds finally an open end such that a corresponding spacetime H exists 
globally and the area dH of its null inner boundary stabilizes at a finite value 47iR B for 
t — > oo. That limiting inner boundary dH is an event horizon and half of Rb is the Bondi 
mass. It follows from the construction that all null geodesies originating in H and directed 
outward are complete. 
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Appendix 

Lemma. 

Assume notation of the main text. Let U, VeH 2 (T,Q Ut ). Define a H k extension Eg"' = 
[xeE : 2m + i] — 5 < |x| < oo] of Eg"' by assuming following initial data in the interval 
[2m + i] — 5, 2m + 77] 

d*f(R) = 0, d k r - l f(R) = d^fiRo), d k r ^}\R) = d^f(R ) - / drd k r ~ i+1 f(r), 

JR 

where f — U,V and i — 2, ...k. Then 

for sufficiently small 5 and some constant C. 

Proof. The two norms in question differ by a sum of terms of the form 
fx°_ s (d k ~ 1 f(R)) 2 R 2 dR; that sum is a polynomial in 5 that can be bounded from above 
by C5yj^ 1 ((9*/(i?o)) 2 - By one of Sobolev inequalities (see the proof of Proposition 7.3), 
\d l r f(Ro)\ < C'\ |i9* +1 /| |sgut; combining all information one arrives at the desired inequality. 
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The initial values of the scalar field in the extension can be obtained from integrating 
U and V, <j>{R) = <p(R ) - Jr° dr±(U + V). That would lead to an estimation of the 
contribution of the potential term to the mass m that is due to the enlargement of the 
initial open end. Together with the preceding Lemma, that would give m(i? — 5) < m for 
sufficiently small 5. 
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